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I discuss an approach to derive from first principles, a real-time formalism to 
study the dynamical interplay of quantum and statistical-kinetic properties of 
non-equilibrium multi-parton systems produced in high-energy QCD processes. 
The ultimate goal (from which one is still far away) is to have a practically appli- 
cable description of the space-time evolution of a general initial system of gluons 
and quarks, characterized by some large energy or momentum scale, that expands, 
diffuses and dissipates according to the self- and mutual-interactions, and even- 
tually converts dynamically into final state hadrons. For example, the evolution 
of parton showers in the mechanism of parton-hadron conversion in high-energy 
hadronic collisions, or, the description of formation, evolution and freezeout of a 
quark-gluon plasma, in ultra-relativistic heavy-ion collisions. 



1 Introduction 

In general, the study of a high-energy multi-particle system and its quantum 
dynamics involves three essential aspects: first, the aspect of space-time, ge- 
ometry and the structure of the vacuum; second, the quantum field aspect 
of the particle excitations; and third, the statistical aspect of their interac- 
tions. These three elements are generally interconnected in a non-trivial way 
by their overall dynamical dependence. Therefore, in order to formulate a 
quantum description of the complex non-equilibrium dynamics, one needs to 
find a quantum-statistical and kinetic formulation of QCD that unifies the 
three aspects self-consistenthj. The main tools to achieve this are: the closed- 
time-path (CTP) formalismlHo (for treating initial value problems of irreversible 
systems), and (ii) transport theory based on Wigner function techniques (for 
a kinetic description of inhomogenous non-equilibrium systems). 

The common feature of high-energy particle collisions is that they allow 
a distinction between a short-distance quantum field theoretical scale and a 
larger distance statistical-kinetic scale, which is essentially an effect of ultra- 
relativistic kinematics. This advantagous property facilitates the passage from 
exact QCD field theory of coherent non-abelian gauge fields to an approximate 
quantum kinetic theory of an ensemble of incoherent gluons. When described 
in a reference frame, in which the particles move close to the speed of light, 
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the effects of time dilation and Lorentz contraction separate the intrinsic quan- 
tum motion of the individual particles from the statistical correlations among 
them. On the one hand, the quantum dynamics is determined by the self- 
interactions of the bare quanta, which dresses them up to quasi-particles with 
a substructure of quantum fluctuations. This requires a fully quantum the- 
oretical analysis including renormalization. On the other hand, the kinetic 
dynamics can well be described statistical-mechanically by the motion of the 
quasi-particles that is, by binary interactions between these quasi-particles, and 
by the possible presence of a coherent mean color field that may be induced by 
the collective motion of the partons. Such a distinct description of quantum 
and kinetic dynamics is possible, because the quantum fluctuations are highly 
concentrated around the light cone, occurring at very short distances, and de- 
couple to very good approximation from the kinetic evolution which is dictated 
by comparably large space-time scales. As mentioned, the natural two-scale 
separation is just the consequence of time dilation and Lorentz contraction, 
and is true for any lightcone dominated process. In fact, at asymptotic en- 
ergies the quantum fluctuations are exactly localized on the lightcone, and so 
the decoupling becomes perfect. This observation is the key to formulate a 
quantum kinetic description in terms of particle phase-space densities, involv- 
ing a simultanous specification of momentum space and space-time, because 
at sufficiently high energy, the momentum scale Ap of the individual particles' 
quantum fluctutions and the scale Ar of space-time variations of the system 
of particles satisfiy ApAr ^> 1, consistent with the uncertainty principle. 

In what follows, I am guided by the recent paper@ and the related literature 
discussed therein, plus on preliminary results of work in progresstj. For purpose 
of lucidity, I will henceforth confine myself to pure Yang-Mills theory, i.e. 
consider gluons only and ignore the quark degrees of freedom. The latter are 
straightforward to include. 



2 Non-equilibrium techniques for QCD 

2.1 Basics of the closed-time-path formalism 

As proclaimed, the goal is to describe the time evolution of a non-equilibrium 
quantum system consisting of an initial ensemble of high-energy gluons at 
starting time to. In this context, the starting point of non-equilibrium field 
theory is to write down the CTP in-in amplitude Zp for the evolution of 
the initial quantum state \in) forward in time into the remote future, in the 
presence of a medium which described by the density matrix. The amplitude 
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Zp is formally given byl3: 

Z P [J,p] = (in\in)j = Tr{uHt ,t)U(t,t )p(to)}j > W 

where ^7 = (^7 + ,^7~) is an external source with components on the + and 

— time branch. p(to) denotes the The initial state density matrix is denoted 
p(to), U is the time evolution operator, and T (T>) denotes the time (anti-time) 
ordering. Within the CTP formalism the amplitude Zp can be evaluated by 
time integration over the closed-time-path P in the complex t-plane. This 
closed-time path extends from t = to to t — t^ in the remote future along the 
positive (+) branch and back to t = to along the negative (— ) branch, where 
any point on the + branch is understood at an earlier instant than any point 
on the — branch. With Zp defined on this closed-time-path, one may then, as 
in standard field theory, derive from it the Green functions and their equations 
of motion. The differences between the CTP and the standard field theory, 
which are briefly summarized below, arise then solely from the different time 
contour. 

The interpretation of this formal apparatus for the evolution along the 
closed-time path P is rather simple: If the initial state is the vacuum itself, 
that is, the absence of a medium generated by other particles, then the density 
matrix p is diagonal and in ([!]), one has \in) — ► |0). In this case the evolution 
along the + branch is identical to the anti-time ordered evolution along the 

— branch (modulo an irrelevant phase), and space-time points on different 
branches cannot cross-talk. In the presence of a medium however, the density 
matrix contains off-diagonal elements, and there are statistical correlations 
between the quantum system and the medium particles (e.g. scatterings) that 
lead to correlations between space-time points on the + branch with space- 
time points on the — branch. Hence, when addressing the evolution of a multi- 
particle system, both the deterministic self-interaction of the quanta, i.e. the 
time (anti-time) ordered evolution along the + (— ) branch, and the statistical 
mutual interaction with each other, i.e. the cross-talk between + and — branch, 
must be included in a self-constistent manner. The CTP method achieves 
this through the time integtation along the contour P. Although for physical 
observables the time values are on the + branch, both + and — branches will 
come into play at intermediate steps in a self-consistent calculation. 

The convenient feature of this Green function formalism on the closed- 
time path is that it is formally completely analogous to standard quantum 
field theory, except for the fact that the fields have contributions from both 
time branches, and the path- integral representation of the in-in amplitude (Q), 
contains as usual the classical action I [A and source terms J o A, but now for 
both time branches, 



3 



Zp[J+,J-,p\ = J VA + VA~ exp [i ( I[A + ] + J + ■ A+^j (2) 

-i(l*[A-} + J~-A~)] M\p] . 

From this path-integral representation one obtains the n-point Green functions 
G^(xi, . . . ,x n ), which are now 

G a ^- a "(x u ...,x n ) , ai = ±, 

depending on whether the space-time points Xi lie on the + or — time branch, 
and it is possible to construct a perturbative expansion of the non-equilibrium 
Green functions in terms of modified Feynman rules (as compared to standard 
field theory), 

(i) The number of elementary vertices is doubled, because each propagator 

line of a Feynman diagram can be cither of the four components of the 
Green functions. The interaction vertices in which all the fields are on 
the + branch are the usual ones, while the vertices in which the fields are 
on the — branch have the opposite sign. On the other hand, combinatoric 
factors, rules for loop integrals, etc., remain exactly the same as in usual 
field theory. 

(ii) All local 1-point functions, such as the gauge-field or the color current, 

are 'vectors' with 2 components, 

A{x) = = ( j- ) (3) 

Similarly, all 2-point functions, as the gluon propagator and the 
polarization tensor n^„, are 2x2 matrices with components 



^ ( A++ A+ ~ \ tt/ \ ( n++ if 
a(xi,.t 2 ) = ^ A _ + A __ j n( Xl ,x 2 ) = n _+ n _ 



(4) 



Explicitcly, the components of the propagator are 

A^(x,y) = A++(x,y) = -i(TA+(x)At(y)) 

A< v (x,y) = A+-{x,y) = -i{At{y)A-{x)) 

A>„(x,y) = A-+(x,y) = -i(A-(x)At(y)) 

Aj v {x,y) = A~~(x,y) = —i(TA~(x)A~(y)}, (5) 
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where A F is the usual time-ordered Feynman propagator, A F is the cor- 
responding anti-time-ordered propagator, and A > (A < ) is the unordered 
correlation function for xq > yo ( x o < Vo)- I n compact notation, 

^(x,y) = -i{T P A(x)A{y)(y)) , (6) 

where the generalized time-ordering operator Tp is defined as 

T P A{x)B{y) := 9 P (x ,y a ) A(x)B(y) + 6 P (y , x ) B(y)A(x) , (7) 

with the #p-function defined as 

n , n f 1 if xq succeeds yo on the contour P , n . 
up(xo,yo) = < . , ,, _ . (8) 

[0 if xq precedes yo on the contour P 

Higher order products A(x)B(y)C(z) . . . are ordered analogously. Fi- 
nally, for later use, let me also introduce the generalized 5p-function 
defined on the closed-time path P: 

!+<5 4 (x — y) if xq and y from positive branch 
— 5 4 (x — y) if Xq and yo from negative branch 
otherwise 

(9) 

Henceforth I will not explicitly label the +, — components, unless it is 
necessary. Instead a compressed notation is used, in which it is understood 
that, e.g., 1-point functions such as A(x) or J{x), 2-point functions such as 
A M (ir, y) or Tl^^x, y), receive contributions from both + and — time branches. 

2.2 The generating junctional for the non- equilibrium Green functions 

The amplitude Zp introduced in ([!]) admits a path-integral representation 
which gives the generating functional for the OTP Green functions defined on 
closed-time-path PEl: 

Zp[J,p] = M J EUdetJF 6{f[A}) exp {* ( / [A, J]) } M(p) , (10) 

where A — (A + ,A~) and J — (J + , J~) have two components, living on the 
+ and — time branches. 

The structure of the functional Zp in (fL0|) is the following: 
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(i) The functional integral (with normalization TV) is over all gauge field con- 
figurations with measure VA = Yl a '^ >y tfi, subject to the condition of 
gauge fixing, here for the class of non-covariant gauges defined by 

f a [A] := h-A a {x) - B a (x) =^ (UPA^ix)) = 0, (11) 

where n M = and n M is a constant 4-vector, being either space-like 

(n 2 < 0), time-like (n 2 > 0), or light- like (n 2 = 0). With this choice of 
gauge class the local gauge constraint on the fields A®(x) in the path- 
integral ( |To| ) becomes, 

dctF 5 (h ■ A a - B a ) = const x exp j-^- J d 4 x [h ■ A a {x)} 2 ^ 

= I GF [h-A] , (12) 

where det.F is the Fadeev-Popov determinant (which in the case of the 
non-covariant gauges turns out to be a constant factor), and where 5{n ■ 
A) = Yl a S(n ■ A a ). The right side translates this constraint into a the 
gauge fixing functional Igf- The particular choice of the vector and of 
the real-valued parameter a is dictated by the physics or computational 
convenience, and distinguishes further within the class of non-covariant 



gauges 



homogenous axial gauge : n 2 < a = 

inhomogenous axial gauge : n 2 < a = 1 

temporal axial gauge : n 2 > a = 

lightcone gauge : n 2 — a = . (13) 



(ii) The exponential I is the effective classical action with respect to both 
the + and the - time contour, I [A, J] = I[A + ,J + ] - I* [A~,J~], 
including the usual Yang-Mills action I YM = J (1 4 x£ym, plus the source 
J coupled to the gauge field A: 

I [A J] = -\J^xT^(x)T^ a {x) + J p d 4 xJ«(x)A^ a (x) 

= Iym[A] + Jo A, (14) 

where T« v = d£A a v - fl*Aj + g f abc A\A%. 
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(iii) The form of the initial state at t = to as described by the density matrix 
p is embodied in the function M. (p) which is the density-matrix element 
of the gauge fields at initial time to , 

M(p) = (A + (t )\p\A-(to)) = exp(iK[A}) , (15) 

where A^ refers to the + and — time branch at to, respectively. The func- 
tional K, may be expanded in a series of non-local kernels corresponding 
to multi-point correlations concentrated at t = to, 

JC[A] = /C (0) + J d 4 xlC^\ l a (x)A fl - a (x) 

+ lJ p d i xdS^%l\x,y)A^ a (x)A^ b (y)... (16) 

Clearly, the sequence of kernels KS n ^ contains as much information as the 
original density matrix. In the special case that p is diagonal, the kernels 
= o for all n, and the usual 'vacuum field theory' is recovered. 



The path-integral representation (|l(j) can be rewritten in a form more 
convenient for the following: First, the gauge-fixing functional IoF[fi ■ A] is 
implemented by using ([l2|). Second, the series representation ( |l6| ) is inserted 
into the initial state functional M.(p). Third, K.^ is absorbed in the overall 
normalization J\f of Zp (henceforth set to unity), and the external source J in 
the 1-point kernel K.^: 

/C (0) := i lnW , /C (1) := /C (1) + J. (17) 

Then ( |i~0| ) becomes, 

Zp[J,p] => Z P [K] = J VA exp{i(l[A^])} , (18) 

where, instead of (p"4|), 

/ [A, K] = I YM [A] + I GF [n-A}+ K {1) oA + i JC (2) o(AA) + ^ IC {3) o(AAA) + ... . 

(19) 
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3 Separating soft and hard dynamics and the equations of motion 

The first step in the strategy is a separation of soft and hard physics in the 
path-integral formalism with Green functions of both the soft and hard quanta 
in the presence of the soft classical field is induced by and feeding back to the 
quantum dynamics. The basic idea to split up the gauge field Ay appearing 
in the classical action Iym [A] into a soft (long-range) part Ay, and a hard 
(short-range) quantum field a M : 

■W = J-0fe +ik - X A a y(k)6(»-k°) (20) 

+ / (2^ e+ *^ (fc) 6{k °-^ = A > ] + < {X) - 

This is the formal definition of the terms 'soft' and 'hard'. The soft and hard 
physics are separated by a (at this point arbitrary) space-time scale A = so 
that one may associate the soft field Ay being responsible for long range color 
collective effects, and the hard field ay embodying the short-range quantum 
dynamics. Consequently the field strength tensor receives a soft, a hard part, 
a mixed contribution, 

J* v (x) = ( F yA A ] + f a yM + <l>UAa]) (x) . (21) 

Now comes physics input. Consider the following physics scenario: The 
initial state is a (dilute) ensemble of hard gluons of very small spatial extent 
■C A, corresponding to transverse momenta k\ 3> /i 2 . By definition of A, or [i, 
the short-range character of these quantum fluctuations implies that the expec- 
tation value (ay) vanishes at all times. However, the long-range correlations 
of the eventually populated soft modes with very small momenta k\ <C n 2 
may lead to a collective mean field with non- vanishing (A). Accordingly the 
following condition on the expectation values of the fields is imposed: 

( A-{x) > { > o for* S S ( < (X) ] ^ ° f ° r aU< ' (22) 

Furthermore, for simplicity the quantum fluctuations of the soft field are ig- 
nored, assuming any multi-point correlations of soft fields to be small, 

(A%( Xl ) ...^">„))« (A%( Xl )) ... <A£ (*„)>, (23) 

i.e. take Ay as a non-propagating and non-fluctuating, classical field. 
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When quantizing this decomposed theory by writing down the appropriate 
m-m-amplitude Zp, one must be consistent with the gauge held decomposition 
( f2p| ) into soft and hard components and with the classical character of the 
former. = 0, M^J > 0. That is, I restrict in the following to a class of 

non-equilibrium initial states of Gaussian form (i.e. quadratic in the holds) 
and do not consider possible linear force terms. 

Substituting the soft-hard mode decomposition j20| ) with the condition 
( p2| ) into (|lj), the functional integral of the in-in amplitude ( |l8|) becomes: 

Z P [fC] = J VAVa exp{i(l[A] + I [a] + I [A, a])} , (24) 

with a soft, hard, and mixed contribution, respectively^. 

Introducing the connected generating functional for the connected Green 
functions, denoted by Q^ n 

W P [K] = -i In Z P [K] , (25) 

from which one obtains the connected Green functions by functional dif- 
ferentiation, in terms of mixed products of a M and A^ helds 



(26) 

/c=o 



where the superscript (c) indicates the 'connected parts'. Specifically, one finds 

Q (2) ^(x lV ) = (a^(x)a b v (y))P = iZf v {x,y) . (27) 

These relations define the soft mean field A and the hard propagators A. 

The equations of motions for A and A follow now as in usual field theory 
by functional differentiation of the effective action, 



r, 



A, A 



= W P [K] ~ K [1) oA - i/C (2) o (iA + A 2 ) . (28) 



Note that the main approximation at this point is the truncation of the infinite 
hierarchy of equations for the rt-point Green functions of the excact theory, to 
the 1-point function (the soft mean field A(x)) and the 2-point function (the 
hard propagator A(x,y)), with all higher-point functions being combinations 
of these and connected by the 3-gluon and 4-gluon vertices. 
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3. 1 Yang-Mills equation for the 'soft ' mean field 

The equation of motion for the soft field AAx), is given by STp/SA — —JC^ — 
K.( 2 ) oA, from which one obtains, upon taking into account the initial condition 
/C^ 1 ) = 0, the Yang-Mills equation for A: 



-A, ab -b 



Xf.i 



(x) = -'#(*) - j p <?yK {2 \f{x,y)A">"{y), (29) 



where [-D,F] = D F — F D with the covariant derivative defined as D = 

D X (A] = - igA\ and F^ = F^[A] = [~D X , D„] /(-ig). The left hand 
side of may be written as 



—A, ab^ -b 



r \,b, 



(30) 
(31) 



where, upon taking into account the gauge constraint (|1 1|), the —n^h\A does 
not contribute, because = (n • A) = n v A U: and where 



-(1)^1 



"(2) M 



(X) 



(32) 



I / n^^^A^^i,^)^ 6 ^)^-^^) (33) 

J P n 1 



,%3 ) 



x A (xi)A (x 2 )A (x 3 ) 



(34) 



On the right hand side of (|29|) , the current j is the induced current due to the 
'hard' quantum dynamics in the presence of the 'soft' field A: 



J» = ??!)„(*) + 1(2) M + 



(35) 



^ / V^^) A vA,6c (a;i,a;2) 



p.- 



i( 2 )» = -y / ft^i W ( o)°A c t(^ ^1,^,^3) 



i=l 
3 



(36) 
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x A v '\ Xl ) A x °> cd (x 2 ,x 3 ) (37) 

x a^^^a^^.j/s) v {0) ffi x ,Ayi,y2,ys) ■ (38) 

Finally, the second term on the right side of (^9|) is the initial state contribution 
to the current, which vanishes for t = x° > to- 

Notice that the function S on the left hand side of (^9|) contains the non- 
linear self-coupling of the soft field A alone, whereas the induced current j on 
the right hand side is determined by the hard propagator A, thereby generating 
the soft field. 



3.2 Dyson-Schwinger equation for the 'hard' Green function 

The equation of motionfor the 'hard' propagator, A^(x,y), is STp/SA 
IC^ /{2i), from which one finds after incorporating the initial condition K,^ 
0, the Dyson-Schwinger equation for A: 



(a^)" 1 - (a^V)" 1 + rC + n; 



b 



(x,y) = IC^{x,y), (39) 



where A = Arj, is the fully dressed propagator of the 'hard' quantum fluctua- 
tions in the presence of the 'soft' mean field, whereas A( ) is the free propagator. 
The polarization tensor II has been decomposed in two parts, a mean-field part, 
and a quantum fluctuation part. The mean-field polarization tensor II incor- 
porates the local interaction between the 'hard' quanta and the 'soft' mean 
field, 

H-lv{x,y) = n (1) a Ml/ (x,y) + U (2) a ^(x,y) (40) 



n (1 /t(^) = l {sUx,y) J p d i zV$% x (x,y,z)A X - c (z) (41) 

n (2 /t(z,y) = y $>(x,y) J *z<twW$% Xir {x,y,z,w) 

xA X ' c (z)A^ d {w). (42) 

plus terms of order g 3 A which one may safely ignore within the present ap- 
proximation scheme. The fluctuation polarization tensor II contains the quan- 
tum self-interaction among the 'hard' quanta in the presence of A, and is given 
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(2,) (2) 

by the variation of 2-loop part T p , of the effective action, 2iST p J /6A, 

n^(x, y ) = (n (1) + n (2) + n (3) + n^V" (x, y ) , (43) 



nri%,{x,y) = -y J d^dSi W {0) a £L(x,y,x 1 ,y 1 ) A x -' cd (y 1 ,x 1 ) (44) 

x A AA '' cc '(a;i,3/i) A™'^'^, y 2 ) V$$ b v (y 2 ,Vi,v) (45) 
z=l 

x A AA ' 1 cc ' (xi , j/i ) A* 7 * 7 ' ' dd ' (x 2 , j/2 ) 

x A TT '' ee ' (ass, 2/ 3 ) WZZiKfa, V2, yi,y) (46) 

x A^'(x 2 ,z 2 ) A^"> e -f"(z 3 ,z 3 ) V}%'/(zs,Z2,Zi) 
xAo x '>f c '(z uyi )A x °'< cd '(x u y 2 ) V^:( yi ,y 2 ,y). (47) 



Note that the usual Dyson-Schwinger equation in vacuum is contained in (39) 
-(|47|) as the special case when the mean field vanishes, A(x) — 0, and initial 
state correlations are absent, K.( 2 '(x, y) = 0. In this case, the propagator 
becomes the usual vacuum propagator, since the mean-field contribution II is 
identically zero, and the quantum part II reduces to the vacuum contribution. 

4 Transition to qantum kinetics 



The equations of motion ([29|) and (|39j) are non-linear integro-differential equa- 
tions and clearly not solvable in all their generality. To make progress, one 
must be more specific and employ now the details of the proclaimed physics 
scenario, described above. 



4-1 Quantum and kinetic space-time regimes 

The key assumption is the separability of hard and soft dynamics in terms of 
the space-time scale r(/i) oc ~ 1 fm. This implies that one may charac- 
terize the dynamical evolution of the gluon system by a short-range quantum 
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scale r qua -C r(p), and a comparably long-range kinetic scale r^in XL ?" (/■*)• Low- 
momentum collective excitations that may develop at the particular momen- 
tum scale g/j, are thus well seperated from the typical hard gluon momenta of 
the order /u, if g -C 1. Therefore, collectivity can arise, because the wavelength 
of the soft oscillations ~ is much larger than the typical extention of the 
hard quantum fluctuations ~ Notice that this separation of scales is not 
an academic construction, but rather is a general property of quantum field 
theory. A simple example is a freely propagating electron: In this case, the 
quantum scale is given its the Compton wavelength ~ 1/ m e in the restframe of 
the charge, and measures the size of the radiative vacuum polarization cloud 
around the bare charge. The kinetic scale, on the other hand, is determined by 
the mean-free-path of the charge, which is infinite in vacuum, and in medium is 
inversely proportional to the local density times the interaction cross-section, 
~ l/(n g (Ti n t). Adopting this notion to the present case of gluon dynamics, let 
me define r qua and rkin as follows: 

quantum scale r qua : Measures the spatial extension of quantum fluctua- 
tions associated with virtual and real radiative emission and re-absorption 
off a given hard gluon, described by the hard propagator A. It can thus 
be interpreted as its Compton wavelength, corresponding to the typical 
transverse extension of the fluctuations and thus inversely proportional 
to the average transverse momentum, 

Tqua = A ~ — ^— , ( k± ) > fl , (48) 

where the second relation is imposed by means of the definition ( po| ) of 
hard and soft modes. Note that Ac is a space-time dependent quantity, 
because the magnitude of (k±) is determined by both the radiative self- 
interactions of the hard gluons and ther interactions with the soft field. 

kinetic scale rki n - Measures the range of the long- wavelength correlations, 
described by the soft mean-field A, and may be parametrized in terms 
of the average momentum of soft modes (q), such that 

r km = A ~ j^— , ( q± ) < g fi , (49) 

where A may vary from one space-time point to another, because the 
population of soft modes A(q) is determined locally by the hard current 
j with dominant contribution from gluons with transverse momentum 
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The above classification of quantum- (kinetic-) scales specifies in space- 
time the relevant regime for the hard (soft) dynamics, and the separability of 
the two scales r qua and Tkin imposes the following condition on the relation 
between space-time and momentum: 

A < A, or ( k x ) w fi > gn fa ( q x ) . (50) 

The physical interpretation of ( |50| ) is simple: At short distances r qua <C l/(<7/x) 
a hard gluon can be considered as an incoherent quantum which emits and 
partly reabsorbs other hard gluons corresponding to the combination of real 
brcmstrahlung and virtual radiative fluctuatiuons. Only a hard probe with a 
short wavelength < r qua can resolve this quantum dynamics. On the other 
hand, at larger distances rkin pa l/(g/j,), a gluon appears as a coherent quasi- 
particle, that is, as an extended object with a changing transverse size corre- 
sponding to the extent of its intrinic quantum fluctuations. This dynamical 
substructure is however not resolvable by long-wavelength modes > r^m of the 
soft field A. 

Accordingly, one may classify the quantum and kinetic regimes, respec- 
tively, by associating with two distinct space-time points x M and the fol- 
lowing characteristic scales: 

s" = - y» ~ A = — , djt = \ - fi£) ~ gii 

g[i 2 

r» ee \ {x» + y") ~ A = - , fl? = 9£ +d% ~ ix .(51) 

The kinetic scale is therefore g 2 fi 2 : The effect of the soft field modes of A 
on the hard quanta involves the coupling gA to the hard propapgator and is 
of the order of the soft wavelength A = l/(<?/z), so that one may characterize 
the soft field strength by 

9Mr) ~ g/i, gF^(r) ~ <?V , (52) 

plus corrections of order g 2 fj? and g 3 ^ 3 , respectively, which are assumed to be 
small. 

The quantum scale on the other hand is /i 2 , because 

K-l ~ fci > ii 2 » <?V ~ 5 iV , (53) 

and one expects that that the short-distance fluctuations corresponding to 
emission and reabsorption of gluons with momenta k± > /i, are little affected 
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by the long-range, soft mean field, because the color force ~ gF acting on 
a gluon with momentum k± ~ (j, produces only a very small change in its 
momentum. 

Concerning the Yang-Mills equation ( p9| ) , one finds then immediately from 
the above scale relations that both the derivative terms d 2 A and the self- 
coupling terms 5 are of the same order and need to be included consistently in 
order to preserve the gauge symmetry when performing a perturbative analysis. 
Of course, if the field is weak, F ^„ <C 3/x 2 , the nonlinear effects contained in 
the function E of ( p9| ) are be subdominant, so that in leading order of g, the 
color fields would then behave like abelian fields. 

4-2 The kinetic approximation 

The realization of the two space-time scales, short-distance quantum and quasi- 
classical kinetic, allows to reformulate the quantum field theoretical problem 
as a relativistic many-body problem within kinetic theory. The key element 
is to establish the connection between the preceding quantum-theoretical de- 
scription in terms of Green functions And a probabilistic kinetic description in 
terms of of so-called Wigner functionsEl. Whereas the 2-point functions, such as 
the propagator or the polarization tensor, depend on two separate space-time 
points x and y, their Wigner transform utilizes a mixed space-time/momcntum 
representation, which is particularly convenient for implementing the assump- 
tion of well separated quantum and kinetic scales, i.e., that the long- wavelength 
field A is slowly varying in space-time on the scale of short-range quantum 
fluctuations. Moreover, the trace of the Wigner transformed propagaor is the 
quantum analogue of the single particle phase-space distribution of gluons, 
and therefore provides the basic quantity to make the connection with kinetic 
theory of multi-particle dynamics. 

In terms of the center-of-mass coordinate, r — + y), and relative 
coordinate s — x—y, of two space-time points x and y, eq. (|5l|), one can express 
any 2-point function Q(x,y), such as A, IT, in terms of these coordinates, 

5, r-~) = G$(r,8), . (54) 

The Wigner transform Q (r, k) is then defined as the Fourier transform with 
respect to the relative coordinate s, being the canonical conjugate to the mo- 
mentum k. In general, the necessary preservation of local gauge symmetry 
leads to additional constraint, but for the specific choice of gauge jjl]), the 



15 



Wigncr transform is simply 



G(r,s) = J (0ie-* fc - G(r,k) , G(r,k) = J d 4 S e ifc ' s G(r,s) . 

(55) 

The Wigner representation (|5^) will facilitate a systematic identification of 
the dominant contributions of the soft field A to the hard propagator A: First 
one expands both A and Ap in terms of the long-range variation with the 
kinetic scale r (gradient expansion), then one makes an additional expansion 
in powers of the soft field A and of the induced perturbations Ap, ~ flA^. 
On this basis, one isolates and keep consistently terms up to order g 2 /j, 2 A^. 

To proceed, recall that the coordinate r M describes the kinetic space-time 
dependence 0(Arfei n ), whereas s measures the quantum space-time distance 
0{Ar qua ). In translational invariant situations, e.g. in vacuum or thermal 
equilibrium, W(r, s) is independent of and sharply peaked about s M = 0. 
Here the range of the variation is fixed by A = eq. (jiq), corresponding 
to the confinement length « 0.3 fm in the case of vacuum, or to the thermal 
wavelength 1/T in equilibrium. On the other hand, in the presence of a 
slowly varying soft field A with a wavelength A = l/(g/j,), eq. (f49|), the s M 
dependence is little affected, while the acquired r M dependence will have a 
long-wavelength variation. This suggests therefore to neglect the derivatives 
of Q(r, k) with respect to r M of order g/u,, relative to those with respect to 
of order /i. 

Hence one can perform an expansion of the soft field and the hard propa- 
gator and polarization tensor in terms of gradients, and keep only terms up to 
first order, i.e., 

Mx) = A^ (r + |) - Mr) + | • d r Mr) 
MV) = \ (r ~ * Mr) ~ i ■ drMr) 



27 ^ w 2 
A( ) y.v (x, y) = A (O )^(0,s) 

A fJlV {x,y) = A^(r,s) ~ A M „ (0, s) + s • d r A^(r,s) 
Ufj, v (x,x) = H M „(r) + -• drll^u^) 

Il fiv (x,y) = Up V (r, s) ~ 11^(0, s) + s ■ <9 r n M „ (r, s) , (56) 

and furthermore, in order to isolate the leading effects of the soft mean field A 
on the hard quantum propagator A, one separates the mean field contribution 
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from the quantum contribution by writing 

A(r,k) = A m (r,k) = A@(fc) + SA^k) 



(57) 



with a translation-invariant vacuum quantum contribution and a r-dependent 
mean field part, respectively, 



A 



A' 



A=0 



A=0 



a x -a^ = -n 

(58) 

where A( ) denotes the free propagator, and the A the mean-field proagator, 
that is, the free propagator in the presence of the mean field, but in the absence 
of quantum fluctuations. 

Given the ansatz (57), with the feedback of the induced soft field to the 
hard propagator being contained in ^A^ , one can expand the latter in powers 

of the soft field coupling gA, and anticipate that it is at most g times the 
vacuum piece Ara, that is, 



6A m (r,k) = J2 lH r )'4)"%]« * gA(r) ■ d k A m (k) , (59) 



and, to the same order of approximation, dfiSArg, ^ v (r, k) ~ g(d£A )d£A^ flt/ . 
Inserting now into eqs. (p^) and (p9|) the decomposition d57l) with the 



approximation (p9[), and keeping consistently all terms up to order g 2 fi 2 A^, 
one arrives (after quite some journey i) at a set of equations that can be 
compactly expressed in terms of the kinetic momentum rather than the 
canonincal momentum k^ (as always in the context of interactions with a 
gauge field). For the class of gauges gauge (Oh amounts to the replacements 



K M = k M - gA^(r) 



&' 



D, 



% - gdlA (r)d k v 



(60) 

and, within the present approximation scheme, one has K 2 A 3> D r A. The 
result of this procedure is: 



n x ' ab -p b 



(r) = - if.(0 



- -9 j 0f^{-K^ v ir,K) +A^(r,K)K v } (61) 
{K\A^}(K) = d^(K) + \ {ti^ a (K),A™(K)} (62) 
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^ Ur > [A] 



(r,K) = - - [U^(r,K),A^(K) 



(63) 



One sees that the original Dyson-Schwinger equation reduces in the kinetic 
approximation to a coupled set of algebraic equations. Recall that (62) and 
( |63| ) are still 2x2 matrix equations mix the four different components of 
A = (A F , A>, A<, A ¥ ) and of ft = (n i? ,n > ,n < ,n F ). For the following it is 
more convenient to employ instead an equivalent set of independent functions, 
namely, the retarded and advanced functions A ret , A adv , plus the correlation 
function A cor , and analogously II. This latter set is more directly connected 
with physical, jpbservable quantities, and is commonly referred to as physical 
representationu: 



A 



= A 1 



A ; 



A adv = A * _ A < 



A cor = A > + A < 

(64) 

Similarly, for the polarization tensor the retarded, advanced and correlation 
functions are defined as (note the subtle difference to (p5)): 



IT 



n 1 



IT 



IT 



adv 



n 1 



IP 



IT 



-ip 



n< 

(65) 

Loosely speaking, the retarded and advanced functions characterize the intrin- 
sic quantum nature of a 'dressed' gluon, describing its substructural state of 
emitted and reabsorbed gluons, whereas the correlation function describes the 

Sinetic correlations among different such 'dressed' gluons. The great advantage 
of this physical representation is that in general the dependence on the phase- 
space occupation of gluon states (the local density) is essentially carried by the 
correlation functions A > A A K , whereas the dependence of the retarded and ad- 
vanced functions, A re *, A adv , on the local density is weak. More precisely, the 
retarded and advanced propagators and the imaginary parts of the self-energies 
embody the renormalization effects and dissipative quantum dynamics that is 
associated with short-distance emission and absorption of quantum fluctua- 
tions, whereas the correlation function contains both the effect of interactions 
with the soft mean field and of statistical binary scatterings among the hard 
gluons. In going over to the physical representation, one arrives at the set of 
master equations: The Yang-Mills equation (|6l]) reads 

9 /t^^{ + Ap,;^)} (66) 
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and the renormalization ( |62| ) and transport equations (|63|), become □ 



K D r , A 



[A] 



^ 

+ 2 



ReA 



, A^_ 



" i{ nC ° r ' Im ^} 

It A c - r \ 

4 1' M J- ' 



(68) 



where II = II + II, and the real and imaginary components of the polarization 
tensor are denoted by 



M 



2 _ 
[if — 



1 



RelV = - (IT 



jadv\ 



n-*)^ r> ee imiv = i(ir--n- ,,„ 

(69) 

Note also, that are the real and imaginary components of the Hard propagator 
are given by the sum and difference of the retarded and advanced contributions, 
respectively, 



ReA M „ = 



A' 



ado 



ImA 



A ret - A 



ado 



(70) 



The physical significance of the ( |67j ) and (^) is the following: Eq. ( |67j ) deter- 
mines the state of a dressed parton with respect to their virtual fluctuations and 
real emission (absorption) processes, corresponding to the real and imaginary 
parts of the retarded and advanced self-energies. Eq. (|S8|), on the other hand 
characterizes the correlations mong different dressed parton states, and the 
self-energies appear here in two distinct ways. The first two terms on the right 
hand side account for scatterings between quasi-particle states, i.e. dressed 
partons, whereas the last two terms incorporate the renormalization effects 
which result from the fact that the dressed partons between collisions do not 
behave as free particles, but change their dynamical structure due to virtual 
fluctuations, as well as real emission and absorption of quanta. For this reason 
II - * is called radiative self-energy, and n cor is termed collisional self-energy. It 
is well knownEI, that the imaginary parts of the retarded and advanced Green 
functions and self-energies are just the spectral density p — ImA, giving the 
probability for finding an intermediate multi-particle state in the dressed par- 
ton, respectively the decay width T, describing the dissipation of the dressed 
parton. The formal solution of (pTj) and (pq) for the spectral density p is 



P(r, k) 



M 2 



(r/2) 2 



Pm 2 



PT 



(71) 
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describing the particle density in terms of the finite width T and the dynamical 
'mass term' M 2 (which in the 'free-field' case are V = M 2 = 0, corresponding 
to an on-shell, classically stable particle). On the right hand side of (|7l|) , the 
second form exhibits the physical meaning more suggestively in terms of the 
'wavefunction'-renormalization (pj^2 = pr=o) due to virtual fluctuations, and 
the dissipative parts {fir — Pm 2 =o) due to real emission (absorption) processes. 



5 Outlook 

What remains to be done is to solve the set of equations (|6^)-(p8|) which is 
the hardest part. For the case of = = F^ v , this was discussed in Ref. □. 
For the present general case, the coupling between hard gluons (A) and the 
soft field (A), complicates things considerably. A possible iterative scheme of 
solution, which is currently under study may be as follows: 

a) Specify initial condition in terms of a phase-space density of hard gluons at 
time t = tQ. This initial gluon distribution determines Ara(t = to, r, k). 



b) Solve the renomalization equation (62) with A(to,r) = 0, i.e. just as in 



the case of vacuum except that now K = k — gA contains the soft 
field. Substitute the resulting form of AS* and A^j" into the transport 



equation (63) to get the solution for Aj^ r . 



c) Insert A^ r into the right hand side of the Yang-Mills equation (|66|), and 

solve for A. 

d) Repeat from a) but now include the finite contribution from the coupling 

between Ara and A. 
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